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We investigate the multiparticle quantum superposition and the persistence of multipartite entan- 
glement of the quantum superposition generated by the quantum injected high-gain optical para- 
metric amplification of a single photon. The physical configuration based on the optimal universal 
quantum cloning has been adopted to investigate how the entanglement and the quantum coherence 
of the system persists for large values of the nonlinear parametric gain g. 



I. INTRODUCTION 

In recent years, a large number of experiments aimed 
at the verification of fundamental aspects of quantum 
mechanics, such as quantum nonlocahty, have been re- 
ahzed by adopting photon particles mutually interact- 
ing through non-linear optical (NLO) process. In addi- 
tion sophisticated NLO methods have been extended to 
relevant investigations and realizations in the domain of 
the emerging new sciences of Quantum Information (QI) 
and Quantum Communication. In particular one of such 
methods, the quantum injected nonlinear (NL) para- 
metric amplification (QIOPA) of single photon states, 
was particularly fruitful since it was adopted to provide 
the first experimental realization of the quantum cloning 
transformation, a fundamental QI concept This 
one provides the optimal distribution of the information 
contained in a quantum system, e.g. a quantum-bit 
state, or quhit onto a system of higher dimension, M > N 
0,1^. By virtue of the isomorphism existing between any 
logic qubit associated with spin-1/2 and the polarization 
state of a single photon, there it is generally supposed 
that photons, identically prepared in an arbitrary state 
of polarization {\(|)))^ are injected into the amplifier on 
the input mode ki . The amplifier then generates on 
the output cloning mode (C) M > N copies, or clones 
of the input qubit \<j)). Moreover, in the case of mode- 
nondegenerate QIOPA, the device simultaneously gener- 
ates M — N states \4>'^) on the output anticloning mode 
k2 (AC) thus realizing an universal quantum NOT gate 

i 

In the last years, the QIOPA scheme has been at 
the basis of experimental realizations of the 1^2 
universal optimal quantum cloning machine (UOQCM) 
0, 0, S and of the 1^3 phase covariant quan- 

tum cloning machine (PQCM) "^10]. These tests, carried 
out in low power linearized conditions, i.e. with very 
low values of the NL parametric gain parameter g << I, 
were followed recently by a series of OPA works, car- 
ried out in absence and in presence of quantum injection, 
which realized the high gain (HG) spontaneous and stim- 
ulated gene ration of large number of output photons M 
hllhihj. Within this new N —!■ M cloning endeavor, 
a multi-photon superposition entangled state was gener- 
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FIG. 1: (Color on line) Schematic diagram of the single pho- 
ton quantum injected optical parametric amplifier (QIOPA). 
The injection is provided by an external spontaneous para- 
metric down conversion source of polarization entangled pho- 
ton states H|.The losses are simulated by the insertion of a 
beamsplitter (BS) over each propagation mode ki. Inset: 
each beamsplitter couples input modes a\j^, with unpopu- 
lated ones {a'jpf). 

ated, indeed a Schroedinger Cat state El El- By virtue 
of the information-preserving (i.e. coherence-preserving) 
property of the parametric process, this implied the de- 
terministic transferal of the well accessible and easily 
achievable quantum superposition condition affecting any 
input single-particle qubit to a "mesoscopic", i.e. multi- 
particle, amplified quantum-state ll^ . 

In the present article, we investigate theoretically the 
nature of quantum injected optical parametric amplifi- 
cation in high gain regime. In particular we intend to 
investigate the most important property of the process, 
i.e. the entanglement of the output modes in the multi- 
particle condition. We show how this difficult task can 
be undertaken by application of a novel technique, here 
referred to as "pair extraction techniqu^\ adopted to in- 
vestigate the multiphoton states generated by high-gain 
spontaneous parametric down-conversion (SPDC) 

In case of a bipartite entanglement, e.g. established 
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over the output cloning and anticloning modes ki and k2, 
the basic idea is to investigate a two-photon field com- 
ponent "extracted" out of the output multi-photon field 
and then infer the entanglement properties of the original 
field. Conceptually, the basic argument underlying this 
method consists of the impossibility of creating or en- 
hancing the entanglement by any local operation, e.g. in 
this case by induced losses In order to do that, in our 
case the adopted theoretical model takes into account the 
propagation losses leading to the imperfect detection of 
the output field. In the SPDC case the explicit form for 
the two-photon output state has been found to exhibit 
a Werner state structure, i.e., consisting of a weighted 
mixture of a maximally entangled (singlet) state with a 
fully depolarized state T?|. This structure is resilient to 
losses for any value of the nonlinear gain parameter g. A 
similar approach will be applied to investigate theoret- 
ically and experimentally the bipartite entanglement in 
the high-gain QIOPA process. The theoretical model en- 
ables to obtain the explicit form of the two-photon out- 
put density matrix for any value of g. More precisely, 
the application of the method goes as follows. Let's start 
from a polarization entangled pair of photons associated 
with the output modes ki and ky: Figurel. The pho- 
ton created over kr freely propagates and is detected by 
the phototube Dt while the photon associated to mode 
ki is injected into the optical parametric amplifier. As 
a result of the coherence-preserving nature of the am- 
plification process the quantum correlations between the 
initial two photon entangled states are transferred into 
correlations between the photon on mode ky (hereafter 
referred to as the "trigger") and the output field of the 
QIOPA device. 

In order to characterize the field generated over the 
three output modes k^, ki and k2 in high gain regime, 
we explicitly derive the reduced density matrix of the 
three photon states obtained over all the output modes 
adopting the pair extraction technique. Hence we 
shall investigate the characteristics of each bipartite sys- 
tems: trigger-cloning, trigger-anticloning and cloning- 
anticloning modes. Each subsystem has been found to 
exhibit a Werner structure. By adopting this model we 
are able to analyze the general entanglement properties 
of the overall 3 mode system. Furthermore, this leads to 
the demonstration of the persistence of entanglement for 
the two systems trigger-cloning and cloning- anticloning 
for any large value of the gain parameter g, in our opinion 
the main result of the present work 

The article is structured as follows. In Section II we 
present the method to extract from multiphoton states 
generated by amplification process two photons compo- 
nents making use of lossy channel and imperfect detec- 
tion. In Section III we analyze the quantum correla- 
tions between trigger, cloning and anticloning, realizing 
a model for the three qubits density matrix with the same 
projection scheme. Finally, Section IV is devoted to the 
discussions of agreement with the experimental appara- 
tus and the perspectives of the present research. 



II. QIOPA IN THE PAIR-EXTRACTION 
REGIME 

In this Section, we analyze the amplified multiphoton 
quantum injected state with particular emphasis on the 
effect of losses on bipartite correlations. Let us consider 
the experimental scheme reported in Fig.l. A single pho- 
ton (mode ki) is injected into the non- linear (NL) crys- 
tal, typically a BBO (/3-barium-borate), cut for Type II 
phase matching and excited by a sequence of UV (ultra- 
violet) mode-locked laser pulses of wavelength (wl) Ap 
propagating along the mode kp . The relevant modes of 
the NL 3-wave interaction driven by the UV pulses as- 
sociated with mode kp are the two spatial modes with 
wave- vector (wv) k^, i = 1,2, each one supporting the 
two horizontal [H) and vertical (V) polarizations of the 
interacting photons. The QIOPA is A-degenerate, i.e. 
the interacting photons have the same wl's A = 2Ap. The 
injected single photon is provided by an external spon- 
taneous parametric down conversion source of biphoton 
states "i]. 

The Hamiltonian of the parametric down-conversion 
process in the interaction picture reads 

H = iK{a\jjdly - a\yaljj) + h.c. (1) 

where a\- represents the creation operators associated 
to the spatial propagation mode k^, with polarization 
j = {H,V}. The NL coupling constant n depends on 
the nonlinearity of the crystal and is proportional to the 
amplitude of the pump beam. It has been theoretically 
shown 0, 01 , that H is invariant under simultaneous 
general SU(2) transformations of the polarization vec- 
tors for modes ki and k2. We may, then, cast the above 
expression in the following form: 

H = iK (a]J>l,^ - al,J)l,'j + h.c. (2) 

where the field labels refer to two mutually orthogonal 
polarization unit vectors for each mode, ^' and cor- 
responding respectively to the state vectors: j^*) and 
I'l'^) . This Hamiltonian generates a unitary transforma- 
tion U — e^'Ti" acting on the input multimode photon 
state. 

Let's consider first the injection of a polarization en- 
coded qubit in the QIOPA system: 

|vI'),^^a|vl/)f^+/3|vE-)}'^ (3) 

with \af + — 1, defined in the 2 x 2-dimensional 
Hilbert space of polarizations (Ir). We analyze this 
configuration by accounting for the 2 interacting opti- 
cal modes ki and k2, i.e. in terms of the basis vectors: 
I*>f7v = Il)ml0)iy 10)2^10)2^ ^ |1, 0,0,0), |vl/)^^ = 
|0, 1,0,0). In virtue of the general information pre- 
serving property of any nonlinear (NL) transformation 
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of parametric type, the output state is again expressed "reflected modes" {kP*^^}. We assume further that all 
by a " multipartide qubit": BS's are equal and characterized by a common spatially 



I*) 



(4) 

Since |^')/jv ^ pure state and U is unitary, j^*) is also a 
pure state. Furthermore and most important, it consists 
of a quantum superposition of two orthonormal, multi- 
particle states \'^)^ and 1^*)^. This is indeed the dom- 
inant property underlying the celebrated Schroedinger 
Cat concept PITH: 



1*)^ ^ 7 E (-r)*+n-iFVTTTN + i,j,j,z) (5) 



j=o 



7 J2 (-rr+-'(-iFVJ+iN,j + i,j,*)(6) 

i, j=0 

where C = coshg, 7 = C^"^, F = tanhg, and g = ntint 
expresses the NL gain of the parametric process, Unt be- 
ing the interaction time. 

In order to implement the pair extraction technique de- 
scribed above in Section 1, we may now explicitly derive 
the theoretical expression of the pair extracted^ density 
matrix in regime of photon losses on all QIOPA spatial 
and polarization output modes. According to a standard 
procedure, the losses are simulated by the insertion of 
dummy ideal beam splitters (BS) on the optical modes, 
followed by ideal detectors {23|. Let us first express the 
above QIOPA output wavefunction in absence of losses, 
in terms of vacuum states: 



i, j=0 



til 



i^'^^^H^ - ^H^^'^^^ 10,0,0,0) 

(7) 

After insertion of the BS's , the lossless QIOPA 
output field operators |a|j|are transformed into the 

l^ij-oc/T l^'^^i^S again on the output lossy channels, by 
the unitary BS map: 



'hj-OUT (*) 



Here the operators |6|j|and |&ij_o(7T| correspond to 
the input and output side, i.e." reflected modes" of the 
dummy BS's that do not coincide with the QIOPA states. 
Precisely, the set |&|^|and |&ij_o(7j-|act correspond- 
ingly on a set of side vacuum states and on the set of 
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g = l 




(c) 



g = 3 f 
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FIG. 2; ((Color on line) Theoretical density matrix of the 
reduced two-photon output state over the output modes ki 
{i — 1,2} conditioned by the injection of a generic input 
polarization qubit \(f>) (left column) and (right column) 

on the mode ki. (a) corresponds to g = 0.1, (b) to g = 1, (c) 
to g = 3. The imaginary parts with all matrix elements equal 
to are not reported. 



and polarization independent transmittivity parameter 
V- 

The lossy effect of the BS's maps \'^) into an output 
state \'^)ouT which involves the four output "transmit- 
ted modes" and the four output "reflected modes" of the 

B S 

two BS's. Precisely, \'^)out expressed in terms of Fock 
states: \niH , niv , ri2H , n2v) a \niH,niv,n2H,n2v)b 
where the two terms in the tensor product repre- 
sent respectively the output transmitted BS's modes 
(a— modes) and the output reflected modes (6— modes). 
The output state of the overaU {QIOPA + BS's} sys- 
tem, is found to reproduce the quantum superposition 
behavior: 
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I*,-, . , ± hli;;*^ (-.yrrr?)— ' („ 1.5) , , ic,)) (s, 

where 

^ ££££ (^)'"™^Vcr)a)(oa)"-"" 

'^■■'^ ^ S££S ( wfc?) ''/alMOO)^ 



with = —1. AU undetected BS "reflected modes" 
must then be discarded, i.e. traced out in the 
above expressions leading to the reduced density matrix 



P' 



over the transmitted modes 



{k^ } that can be expressed as follows: 

oo oo oo oo 
ki=0 k2=0 k3=Q ki=0 



>| \kl,k2,k3, 
(12) 



Up to now we have considered arbitrary, polarization- 
symmetric channel losses. In the following we make the 
additional assumption of high losses, which will greatly 
simplify the calculations. The explicit expression of p' 
may now be easily obtained in the high loss regime, i.e. 
identified by the relation rjn << 1, rjn being the average 
number of photons transmitted by the BS per mode. In 



virtue of this condition we may drop the terms of the 
sums proportional to ry" for n > 2. This corresponds 
to consider only matrix elements of a representation 
in which any Fock basis state correspond to a photon 
occupation number n < 2. As final step we assume to 
detect one photon on each mode by a standard 

2-photon coincidence technique over the QIOPA output 
^■pdes, kiand k2. By this technique we only inves- 
tigate output states involving 1 photon emitted over 
these two modes. This output condition is expressed 
by a matrix representation of p' involving the basis 
states:{|l, 0, 1, 0) , |1, 0, 0, 1) , |0, 1, 1, 0) , |0, 1, 0, 1)}, 
corresponding to the basis: 

{\H), \H)^ , \H)^ \V)^ , \V), \H)^ , \V)^ The 
explicit expression of the detected "pair extracted" 
density matrix is finally given by the expression: 



'd+|/3|'2t2 
al3*d 
-a/3* 2 




Pa*d 



\aV2{l + d)- 
Iap2-|, 
a/3*2 



lap 2 



-/3a*2 



lal^ d - 



|/3|^2(1 
-aP*d 



d) 





/3a* 2 
-I3a*d 
■ 2t2 + \B 



(13) 



'dj 



with d = {1 + t^) and t = V{1 — rf). Let us consider 
the particular asymptotic case 77 — > 0, i.e., t ~ tanhg. 
Figure 2 refers to different injection states and to dif- 
ferent values of the interaction parameter g. There the 
tomographic patterns are reproduced identically for any 
couple of different orthogonal input states { | </)), | </)-'-)} as 
a consequence of the universality of this process, Eq.(|2J). 
Furthermore, the structure of the 4x4 matrices shows 
the relevant quantum features of the output state. For in- 
stance, the highest peak on the diagonals expressing the 
quantum superposition of the input state shifts from the 



position 



to 



in correspondence 



with the OPA excitation by any set of generic orthogonal 
injection states {|0) , j'/'^)}. 



Let us estimate the entanglement of the pair extracted 
output state H13|) . Owing to the tested universality of the 
amplification process we restrict the present analysis to 



the input qubit Y^] 



IN 



\H), i.e. a = 1, /? = 0. In this 
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case the two-photon density matrix reads: 



3(1 + t2) 



+ 

(2 + t2) 

-1 

\ 












The entanglement measurement, the 
C{p') " 



(14) 

concurrence" is 



3(l+t^) 



1 - |\/TTl2)[2il. We found C(p') > 



for any value of g and C(p') ^ for g — > oo. This re- 
sult complies with the one obtained in the case of the 
1 — > M universal quantum cloning by |22l |. There it 
was shown that, in this process closely parallel to the 
present configuration, the concurrence between a clone 
and an ancilla is always different from zero for any value 
of M and vanishes in the limit M — > oo. 

We now analyze the output state (|14|l over each one of 
the output modes, by further tracing p' . For the mode 
ki we obtain: 



Pki 



1 / 5 + 3<2 
6(l + t2)lv 





l + 3t2 



(15) 



The fidelity between the output state and the input is: 
F{\H) , Pi,) = {H\ p',, \H) = In the limit g ^ oo 

we get F = I while the limit g ^ leads to F = |. In 
the present investigation, dealing with the generation of a 
pair of photons, we retrieved exactly the result obtained 
for the 1^2 optimal cloning process. The average num- 
ber of photon over the mode ki is equal to 3n -I- 1 with 
n = sinh^ g. Hereafter the mode ki will also be referred 
to as the cloning mode. 

The single photon polarization state over mode k2 is 
described by the density matrix for any g value: 



, _ / 1/3 
Pk2-\ 2/3 



(16) 



The fidelity between the output state and the orthogonal 
of the input one F{\V) ,^'^2) — 2/3. Hereafter the mode 
k2 will be referred to as the anti- cloning mode. 



III. SCHROEDINGER CAT STATE: 
PERSISTENCE OF ENTANGLEMENT 

Let's now consider a different, more complex configu- 
ration in which the injected photon belongs to a polariza- 
tion entangled pair: see Fig. 3. The previous theoretical 
method will be here adopted to investigate the proper- 
ties of the overall output field. Let us now go into details. 
The initial pump beam is split by an unbalanced beam- 
splitter (BSp) into two beams. A low intensity beam k'p 
excites the NL crystal A which generates pair of entan- 
gled photons over the modes ki and k^: 



kT,kl 



'H\H)kT\H)ki-\V)kT\y)ki) (17) 

where subscripts ky, ki refer to trigger and injection pho- 
tons, propagating along ky and ki modes, respectively. 




NL 

crystal A T 



FIG. 3: (Color on line) Schematic diagram of the quantum in- 
jected optical parametric amplifier (QIOPA) . The pump beam 
k'p excites the NL crystal A which generates pairs of entan- 
gled photons over the modes ki and kr. The photon over 
the mode kr provides the trigger of the overall experiment 
while the twin photon over the mode ki excites the NL crys- 
tal B together with the pump beam kp 4]. The losses are 
simulated by the insertion of a beamsplitter (BS) over each 
propagation mode ki. 



The pump power is sufficiently low to avoid the simulta- 
neous generation of more than one pair of photons. The 
photon over the mode k^ provides the trigger of the over- 
all experiment while the twin photon over the mode ki 
excites the NL crystal B together with the pump beam 

kp y. 

The overall dynamic of the process is described by the 
unitary operator V = IkT®Uki acting on the initial state 
1^ )Tfei ® |0)fe2' where Uki is the time evolution operator 
acting on the injection state and IkT is the unit matrix 
acting on trigger state. Using the results of the previous 
section we obtain 



1^) - ^\^')kT,kl= (18) 

00 

2-^{l^)feT ® E (-r)^+-'(-i)-''VTTT|z + i,j,M) 

00 

- \y)kT ® E (-r)^+n-i)^ V7TT|.,j + l,J,^)}20) 

This multiparticle quantum state exhibits an entangled 
structure connecting the microscopic property of the sys- 
tem, i.e. the single particle "trigger" acting on the 
mode k^, and the macroscopic quantum superposition. 
This indeed corresponds to the original definition of 
" Schroedinger Cat State" [11 0. 

In the following we apply the pair extraction method to 
analyze the entanglement properties of the overall wave- 
function (|19|l . High losses are introduced over the modes 
ki {i — 1.2}. Detection of photons over the three modes 
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ki , enables to reconstruct the density matrix through 
a tomographic technique. The theoretical model shows 
the presence of entanglement between the trigger pho- 
ton and each one of the amplified photons. Since the 
entanglement is not ascribable to the subsequent ampli- 
fication process which is acting locally on the ki arm 
of the pair, one is forced to conclude that the original 
trigger-injection entanglement has survived the QIOPA 
amplification. 



FIG. 4: (Color on line) Theoretical density matrix of the 
reduced three-photon output state over the modes kr, ki {i — 
1, 2} (a) corresponds to g = 0.1, (b) to g = 1. The imaginary 
parts with all matrix elements equal to are not reported. 



Under the pair extraction approximation on k^ {i ~ 
1.2}, i.e. after applying high losses on the output modes 
before detection, the expression of the overall normalized 
3 qubit pair extracted density matrix is easily found: 
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(21) 



Figure 4 reports the various density matrices for different 
values of interaction parameter g. 

In order to understand at a deeper level the multipar- 
ticle quantum cloning process, we may consider the 4x4 
dimensional matrices obtained by tracing p" over the re- 
maining state of the manifold {kj, ki^} . 

Let us start from the system consisting of the kj 
{i = 1.2} modes, i.e. by the cloning-anticloning sys- 



tems. The reduced density matrix p'^^ — TrkT {p") of 
the extracted 2-photons is found: 



I A 



Pkl.k2 





V 





1+p 







p 

_p 1+p 
"2 4 










(22) 
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with p 



We note that the above density 

Pw 

i.e. 



2 1 

3 T+F 

matrix has the form of a Werner state (WS) 

+ ^^I, i.e. a mixture of the maximally 
entangled state |^'-)i.2 = 2-i/2(|i/)^ \V)^ - \V)^ \H)^) 
with probability p and of the fully chaotic state //4, be- 
ing / the identity operator on the overall Hilbert space. 
Note that for any finite value of the NL value g, the 
singlet probability is p > i reaching asymptotically the 
value p = i for 5 — > oo. Since the condition p > i is 
a necessary and sufficient one for state non-separability 
of any WS , the entanglement condition affecting the 
pair extracted p'^^ ^2 1 ^^'^ then of its multiparticle coun- 
terpart pki,k2, does persist for any finite value of g. All 
this can be compared with a similar result obtained in 
regime of spontaneous parametric down-conversion, i.e., 
with no injection of a single photon state. In the SPDC 
case it was found: p = jq^^ |l3j |. The concurrence of 

the state is Cki,k2 = 2 max ^ ^ ^^j ~ 2 
i(sinh^ g + cos^ g)^^- In case of high gain value we find 
Cki,k2 ^ = and TkT.ki = T^T— = =7- The average num- 
ber of output qubits for each mode is equal to ~ Sn. 

A similar analysis can be carried out for the correlation 
affecting the modes and ki , as shown by the structure 
of the pair extracted PkT ki = Trk2 {p")'- 



PkT,kl 








V 





4 









1-9 
4 










1+9 
4 



(23) 



with q = Iy^t^-. Note that once again the density ma- 
trix bears the WS structure p = q\^-) ($-| + ^-f, 
i.e. is a mixture with probability q of the maximally 
entangled state |$-)i^2 = 2-i/2(|i/)^ \H)2 - \V)^ \V)2) 
and of the maximally chaotic state I /A. Once again 
the entanglement condition (7 > ^ is met for any finite 
value of g and g — *■ i for g ^ 00. Note also that in 
the limit of small g w is: g sa | (the present ap- 
proach considers at least the generation of a pair of pho- 
tons). The concurrence of the state of Ea. (|23|) is found: 
CkT,ki = ^(sinh^g -I- cosh^g)"^. In case of large gain 

value is found: CkT,ki ^ = and TkT,ki = C'fcT.fci = ^ 
The average number of output clones is equal to 
3rT + 1 ~ 3n. 

Consider at last the pair extracted reduced density ma- 
trix involving the k2 and kj- modes: p'kT k2 = ^'"fei {p")- 
Its expression is found to be gain independent 



// 

PkT,k2 



i 

i 

^ i 

6 6 



(24) 



It also exhibits a WS structure, p'lx k2 — 



/ with I — 



\k2 - ' \^+/ K'^+l 
i and |$+)i_2 = 2-i/2(|iJ)Jff)2 



system ki- k2 system ki- kx 




(b) 

g=l 0.2\ 






FIG. 5: ( Figure 5. (Color on line) Theoretical density ma- 
trix of the reduced two-photon output state over the output 
modes ki, k2 (left column) and ki, kr (right column), (a) 
corresponds to g = 0.1, (b) to g — 1, (c) to g = 3. The 
imaginary parts with all matrix elements equal to are not 
reported. 



1^)1 1^)2)- ^ consequence, the above state is sepa- 
rable, a result in agreement with the theory of the op- 
timal quantum machines |2^. Indeed the QIOPA has 
been shown to implement over the mode k2 the uni- 
versal optimal dipping machine 0. The completely 
positive (CP) map which implements the UNOT gate 
has the following Kraus representation £unot{p) = 
i {axpcrx + (JYPCTY + <yzP<yz)- The state Pkj'k2 can be 

the 



lkT,kl 



obtained applying to the entangled state |$ 
identity over the mode ky and the CP map Sijj\iQx over 
the mode ki: 



LkT ' 



^-UNOT-kl )kTM(^ IfcT.fcl) 



lkT,k2 



\kT.k2 



I*- 



kT.k2 



lkT,k2 
\kT,k2 



\kT,k2 ' 



(25) 

Interestingly enough, the last expression can be shown to 
be equal to Eq. I|24|) . This confirms the overall validity of 
the present approach. 
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IV. DISCUSSIONS AND CONCLUSIONS 

In the present paper, theory of quantum injected op- 
tical parametric amplification has been extensively in- 
vestigated in regime of high-gain and high losses with 
particular attention for the entanglement properties of 
the output fields. We have exploited the developed the- 
oretical tool to investigate the properties of entangled 
states after a cloning process, demonstrating the persis- 
tence of entanglement for any clone-trigger subsystem. 
Connections with Werner state have been established in 
the physical process of stimulated emission. In addition, 
we have shown that the properties of QIOPA output 
state do comply with the ones conventionally expected 
for a Schroedinger Cat system 0, 0. With respect 
to the persistence of the coherence of the latter system 
for increasing " size" , an interesting problem could be the 
investigation on the entanglement persistence for increas- 
ing the value of the gain g, i.e. for a very large number 
of generated clones. Unfortunately we expect that, for 
an increasing number of clones, the small amount of bi- 
partite entanglement should not be observable due to ex- 
perimental imperfections (like walk off effects and other 
sources of decoherence) . However the density matrices 
experimentally measured in Ref.flll [T3. have been 
found in good agreement with the theoretical ones. On 
the other hand, the technique introduced above turns out 
to be an useful tool to investigate single photon features 
of mesoscopic fields. 

An interesting aspect which deserves further investiga- 
tion is the effect of losses on the orthogonality condition 
of the two wave functions of Eq.5 and 6, in particular how 
the detection efficiency influences the distinction of the 
two initial orthogonal terms. The two extreme condition 
can easily been derived. While for vanishing losses the 
Hilbert-Schimdt distance of the two interfering state is 



found to be dd^-)^,!*)^) = Tr[{\^)" (^H^)^ {^\)^] = 
2, the condition of very low efficiency of detection leads to 
d{p'^, p'^) = 2/9. Finally a new approach to investigate 
quantum features of optical field based on the combina- 
tion of data obtained with different detection efficiencies 
|2^ could improve the present results. 

The theoretical results obtained above have been care- 
fully tested by the experiment reported in Ref.J^]. The 
adopted optical scheme was similar to the diagram shown 
in Figure 3 but for a more compact "folded" configu- 
rations by which a single NL crystal slab, excited in 
both directions by the UV "pump" laser beam, real- 
ized in sequence the SPDC and the QI-OPA operations. 
The experimental investigation of the multiphoton su- 
perposition and entanglement implied by Eqs [2] 
was carried out by means of quantum state tomography 
(QST) according to the pair extraction method previ- 
ously described. The beams associated with the output 
modes ki(i=l,2) were highly attenuated to the single- 
photon level by the two low transmittivity BS's. Ex- 
perimentally the maximal value of gain obtained has 
been found ^oxp = (1.19 ±0.05), while the detection 
quantum efficiencies: read 77i=(4.9 ± 0.2)%; 7?2 = (4.2 ± 
0.2)%. By the previous values we found the condi- 
tion r/n ~ 0.1. The QST analysis of the reduced out- 
put state p'li i,2 determined by the set of input |^')i„: 
{\H), \V), |±)}. The good agreement between theory 
Pki k2 ^^'^ experiment p'li^jF2 is expressed by the mea- 
sured average Uhlmann "fidehty": ^ {Pki^k2^ p'ki k2) = 

\Tr{^^,p'L,U24^2V^? = {^^-^ ± l-2)%. Finally, 
the QST reconstruction of the density matrix p"'^^P has 
been carried out with a fidelity: J^(p"'=''P, p")= (85.0 ± 
1.1)%. 

Work supported by Ministero dell'Istruzione, 
deU'Universita e della Ricerca (PRIN 2005). 
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